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We theoretically investigate the effect of various symmetries of the superconducting order pa- 
rameter A(lo) on the normalized tunneling conductance of SIN junctions by directly solving 
the real-axis Eliashberg equations (EEs) for a half-filled infinite band, with the simplifying as- 
sumption /i* — 0. We analyze six different symmetries of the order parameter: s, d, s + id, 
s + d, extended s and anisotropic s, by assuming that the spectral function a 2 F(Q,) contains an 
isotropic part a F(Q)i B and an anisotropic one, a F(ft) an , such that a F(Q)a. n =g ■ a F(Q)i s , 
where g is a constant. 

We compare the resulting conductance curves at T = 2 K to those obtained by analytical 
continuation of the imaginary-axis solution of the EEs, and we show that the agreement is not 
equally good for all symmetries. Then, we discuss the effect of non-magnetic impurities on the 
theoretical tunneling conductance curves at T = 4 K for all the symmetries considered. 

Finally, as an example, we apply our calculations to the case of optimally-doped high-T c 
superconductors. Surprisingly, although the possibility of explaining the very complex phe- 
nomenology of HTSC is probably beyond the limits of the Eliashberg theory, the comparison 
of the theoretical curves calculated at T=4 K with the experimental ones obtained in various 
optimally-doped copper-oxides gives fairly good results. 

I. INTRODUCTION 

The semi-phenomenological Migdal-Eliashberg theory JlJ has been successfully used in the past to describe many 
properties of low-T c superconductors. For example, it has provided a quite precise explanation of the tunneling 
experimental data obtained on almost all the conventional superconductors. In most of these cases, theoretical 
predictions in good agreement with the experimental results have been obtained by solving the Eliashberg equations 
(EEs) for an s-wave order parameter. However, other symmetries of the order parameter may exist. It is therefore 
interesting to study how, in the framework of this theory, the gap symmetry influences the theoretical conductance 
curves. 

In this paper, we will calculate the theoretical normalized conductance curves of SIN junctions for various gap 
symmetries (s, d, s + id, s + d, anisotropic s and extended s) [|| by directly solving the real-axis Eliashberg equations 
(EEs) in the half-filling case. 

Incidentally, this procedure is much more complicated than the usual one, which consists in solving the EEs in the 
imaginary-axis formulation and then continuing the solution to the real axis. Nevertheless, it is much more general 
and can be used at every temperature, whilst the second approach is meaningful only at very low temperatures (in 
principle, T — > 0). Moreover, we will show that, at a finite but still very low temperature (T = 2 K) the imaginary-axis 
procedure gives results whose agreement with those directly obtained from the real-axis EEs is not equally good for 
all symmetries. 

We will also present the theoretical conductance curves obtained at various temperatures, and with increasing 
values of the coupling constant. It will be shown that large values of the isotropic coupling constant produce some 
characteristic and well-recognizable features of the normalized conductance curves, which are actually observed ex- 
perimentally. 

Moreover, we will investigate the effect of different amounts of non-magnetic impurities, in both the unitary and 
non-unitary limits, on the tunneling normalized conductance curves, and discuss the results for all the symmetries of 
the order parameter. 

Finally, we will try a comparison of the theoretical tunneling curves with experimental data appeared in literature 
and concerning some high-T c compounds. 
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II. ELIASHBERG EQUATIONS FOR DIFFERENT PAIR SYMMETRIES 



In this section we calculate the theoretical normalized conductance for a tunnel junction within the framework of 
Migdal-Eliashberg theory. To do so, we solve the generalized real-axis EEs jjl ^, [| for the renormalization function 
Z(<jj,li) and the order parameter A(w,k). As well known, in the real-axis formalism the EEs take the form of 
a set of coupled integral equations, whose kernels contain the retarded electron-boson interaction a 2 F(fl, k, k'), the 
Coulomb pseudopotential fi* (k, k'), and the energy of the carriers £k, measured with respect to the bare band chemical 
potential |^, 0). In the following, we will use a single-band approximation, and we will restrict our discussion to the 
2-dimensional case, thus referring, for example, to the ab planes of layered superconductors and neglecting the band 
dispersion and the gap in the c direction. For simplicity, we will assume that the Fermi surface is a circle in the ab 
plane j^] and that the wavevectors k and k' are completely determined by the respective azimuthal angles <fi and <f>' , 
since their length is, as usual, taken equal to kp. 

To allow for different symmetries of the pair state, we expand the superconducting spectral function a 2 F(Q, 0, <p') 
and the Coulomb pseudopotential fi*(4>, <p') in terms of basis functions "0is and -0 an , where the subscripts mean isotropic 
and anisotropic, respectively. At the lowest order we suppose a 2 F(Q, </>, tfi') and i^*(4>, <f>') to contain separated isotropic 
and anisotropic contributions, as in the following expressions: 

a 2 F(Q, 0, </>') = a 2 F(n) is Tjj is (0) ^ (</>') + « 2 F(fi) an ^ an (0) </> an (</>') (1) 

</>') = (0 *I*B (0') + iCVW (4>) ^a„ (</>') • (2) 

The basis function (cf>) and ipan (4>) are chosen as follows: 

V'is fo) = 1 

%/2cos(20) (d-wave) 

V'an {4>) — \ 8-\/2/35cos 4 (2(f)) (anisotropic s) (3) 

2^2/3 cos 2 (2(f)) (extended 



We search for solutions of the EEs having the form: 

A(w, 0) = A is (cj) ± A an (w)^ a „ (4>) (4) 
Z(oj, (j>) = Z is (uj) ± Z an (w)?A an (0) . 

We assume as an "ansatz" that the negative sign in the above expressions is used only for the extended s-wave 
symmetry, while in all other cases the positive sign is chosen. In doing so one recovers in the BCS limit the form 
of A given in Ref. 0]. Notice that the choice of the sign in the expression for Z(ui,<p) has no effect. In fact, using 
cqs.([|) makes the Eliashberg equations for A(u,(f>) and Z(u,(j>) split into four equations for A is (u;), Z- is (lu), A an (w) 
and Z an (uj). The equation for Z ari (uj) is a homogeneous integral equation whose only solution in the weak-coupling 
regime is Z an (uj) = 0. Even though in the strong-coupling limit a non-zero solution could exist above a certain 
coupling strength threshold, we do not consider here this rather exotic case and then we assume that the stable 
solution corresponds to Z &n (oj) = for all couplings [S, II. Consequently, we will not write the equation for Z &n (oj). 

In writing the remaining equations for A; s (o;), Z\ s (lJ), A an (w), we insert an additional term which takes into account 
the presence of non-magnetic impurities jl0| . Actually, in the following numerical solutions we will consider this term 
only where explicitly specified, and we will disregard it otherwise. The EEs are then written in the form: 

A i8 (w)^ s (w) = — / dcf> / BB[P a (u\4>)][K^(u,u\<f>)-^e(u B -u;')]du' 
27T Jo Jo 

+i^r — 2 — =2 — (5) 

c2 + p an M + p is M + A is H 

^ p2ir p+oo 

A an (uj)Z is (uj) = — d0 / Re [F an (w', 0)] [K an+ (uj, uj' , 0) - /x* n (w c - u') tp an (<p)] da/ 
^ Jo Jo 

+i^r — 2 — — 2 — (6) 

c 2 + p an H + p is M + A is H 

1 p2lX p+OO 

[1-Z„(w)]e> = — J d0 jf Re[iV is (u/,0)]# is ^u;,u/,<£)du/ (7) 



c 2 + pIm + pIm + nI(u J ) 
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where u> c , which appears in eqs.(|], 0), is a cutoff energy, and the functions P s (u,</>), P an (w,</>) and Ni s (u>,<f)) are 
defined as follows: 

PisM)= ; AisM (8) 



p / ^ A an (^)V> an (0) 

P an (u,0) = (9) 

^/w 2 - [A iB (w) + A an (^ a „ (0)] 2 
N is (u;, (f>) = = = (10) 



[A is (w) + AanMV'an (</>)] 



2 



while 



and 



1 



2;r 



P iSl anM = — / P is , an (w,0)d0 (11) 



2tt 



W to (w) = J- / iV is (^0)d0. (12) 



2tt 

The functions -fCi S , an ±(w, u/, 0) in the integral over u>' are given by 

r + rxi 

K is . an ±(u},uj' ,(f)) = / dfi a 2 F(0) iSian V'is.an (0) • (13) 
Jo 

/(-o/)+n(Q) f{-cj') + n{Q) f(u;')+n{Q) /(w') + n(fi) 



Q + lu' + lu + i5 + u;' — oj — i8 f2 — a/ + a; + i<5 f2 — a/ — — i<5 

where / (w) and n (f2) are the Fermi and Bose distributions, respectively. 

As anticipated above, the last term on the right-hand side of equations (H]?]) allows for the contribution of scattering 
from impurities. It contains the parameters T and C which are defined through the relations T = tii/N(0)tt 2 and 
C=cot((5o), where n\ is the impurity concentration, iV(0) is the value of the normal DOS at the Fermi energy and 8q 
is the scattering phase shift Q . 

Notice that, in the case we want to obtain as + id solution, we simply need to replace the denominator of P an (w, (f>), 
P; s (w,</>) and Ni s (uj,(f>) with yju 2 — [A 2 (uj) + A^(w)^ (</>)], where the generic subscripts 'is' and 'an' have been 
replaced by 's' and 'd', respectively. 

To simplify the problem of numerically solving these equations, we put a 2 P(f2) an = g-a 2 F(fl)- ls where g is a constant 
H pTlfl . As a consequence, the electron-boson coupling constants for the isotropic-wave channel and the anisotropic- 

wave one, which are given by Ai S =2 J +O ° dSl a 2 P(f2)i s /f2 and A an =(l/7r) / Q 27r d<^> V an (0) Jo + °° a2 ^(^)an/^ respec- 
tively, result to be proportional: A an = g ■ X ls . 

The real-axis generalized Eliashberg equations (|^-0) are numerically solved by using an iterative procedure, which 
stops when the real and imaginary values of the functions Aj s (w), A an (w) and Z- m {uj) at a new iteration differ less 
than 1 • 10~ 3 with respect to the previous ones. This usually happens after a number of iterations between 10 and 15. 

By starting from the real-axis solutions Aj s (w) and A an (w), we then calculate the normalized conductance G n (V) of 
a tunnel junction, which is given by the convolution of the quasiparticle density of states N(u>) = Re [sgn(u) Ai s (u;)] 
with the Fermi distribution function: 

Gn(V) = - ( + duj N{uj + eV) • [/(«) - J{uj + eV)] . (14) 
e J-oo 

Of course, this is the simplest possible choice. We are aware that, when anisotropic superconductors are involved, 
the normalized conductance should have a much more complex form to take into account its dependence on the 
junction's geometry and on some other factors such as, for example, the direction of the incident current and the 
angle a between the normal to the interface and the crystallographic a axis p2[ . Here we restrict ourselves to the 
simplest case of plane interface and normal current and we also put a = 0. 

The main features of the resulting G n (V) curves will be discussed in the following sections. 

As already pointed out in the introduction, the direct numerical solution of the real-axis EEs described so far is a 
quite complicated and time-consuming procedure. However, as well known, the EEs admit an imaginary-axis formu- 
lation, in which integrals are replaced by summations over the integer index of the so-called Matsubara frequencies 
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Voltage (mV) Voltage (mV) 

FIG. 1: Theoretical normalized conduc- FIG. 2: Comparison of the normalized tunneling conductance curves ob- 

tance curves for fi* = at T = 4 K in vari- tained from the real-axis solution of the Eliashberg equations, with those 

ous symmetries of the order parameter: s d, obtained by the analytical continuation of the imaginary-axis solutions, at 

s+id (a) and s+d, extended s, anisotropic s T = 2 K . in al l tne symmetries considered, 
(b). 



uj n . The numerical solution of the EEs in this formalism is therefore a much simpler task, which gives the real discrete 
functions A (iw n ) and Z(iu) n ), that can be analytically continued to the real axis by means of Pade approximants p3| . 
In principle, this analytical continuation is meaningful only at very low temperatures. Therefore, one could wonder 
if, and in what temperature range, this method is able to give results in accordance with those obtained by starting 
from the real-axis EEs. A comparison between the density of states (DOS) obtained by the two techniques was made 
for the s-wave symmetry and for small values of the coupling constant (A < 2) fl3[ |, and a substantial agreement was 
found. Here we compare the normalized tunneling conductance curves calculated at T = 2 K not only in the s-wave 
symmetry for greater values of A, but also in the other symmetries analyzed here. The results will be discussed in the 
next section. 



III. THEORETICAL RESULTS 

The numerical solution of the EEs (j^)-(0) requires the knowledge of the function a 2 F(Q)i a which appears in the 
expressions of i^ani and Ai Sian . Even though we are not referring to a particular material, we can make a choice 
for a 2 F(ri)i s by taking advantage of the fact that its detailed shape has little influence on the final solution, which 

instead depends, as observed elsewhere [|l4| , on the quantity uj\ og = exp J +o ° d£l a F ^f^ In fij . Since in this section 

we are only interested in calculating theoretical curves, we can take: 

a 2 F(n) is = ?j-a 2 F(n) Bi2212 . (15) 

Here a 2 -F(r2)f3i22i2 is the electron-boson spectral function determined in a previous paper by inversion of the s- 
wave EEs starting from tunneling data obtained on Bi 2 Sr 2 CaCu208+5 (Bi-2212) break junctions []l5), and A is the 
corresponding coupling constant. The values of Ai s and that of A an are chosen to have T c = max(T c ls , T* n ) = 95 K 
(which, in most of our calculations, is chosen as a representative value for high-T c superconductors). 

From now on, we will put the Coulomb pseudopotential /j,* to zero in all the equations of the previous section. It 
will be shown in the next section that this does not limit the generality of the results. 

As already pointed out, we search for a solution of the real- axis Eliashberg equations which contains an isotropic 
and an anisotropic part, as shown in equation (|]). Actually, the choice of the coupling constants A; s and A an to 
be used for the numerical solution can affect the symmetry of the order parameter A(cl>,</>), which can show either 
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Voltage (mV) Voltage (mV) Voltage (mV) 

FIG. 3: Theoretical normalized con- FIG. 4: Theoretical normalized con- FIG. 5: Theoretical normalized con- 
ductance curves for fj,* = at T=2, 4, ductance curves for /i* — at T—2, 4, ductance curves for /i* = at T—2, 4, 
40 and 80 K in the s-wave symmetry 40 and 80 K in the s+id-wave symme- 40 and 80 K in the anisotropic s-wave 
(a) and in the d-wave symmetry (b) . try (a) and in the s+d- wave symmetry symmetry (a) and in the extended s- 

(b). wave symmetry (b). 

pure isotropic-w&ve or anisotropic-w&ve symmetry, or a mixed symmetry. Moreover, for some particular values of the 
coupling constants, the final symmetry of the order parameter depends on the starting values of A; s (o;) and A an (u;). 

Results with the expected symmetry, and compatible with the value of the critical temperature we have chosen above 
(T c = 95 K) are obtained with the values of the coupling constants reported in Table I, where also the corresponding 
isotropic and anisotropic critical temperatures are shown. 

It is worthwhile to remark that, in all cases apart from the s-wave one, different couples of A; s and A an values can 
give rise to the same critical temperature and the same symmetry and, therefore, Table I only presents one of the 
possible choices. These coupling constants will be used in all the calculations presented in this section. 

Figure 1 reports the theoretical normalized tunneling conductance curves G n (V) at T = 4 K in the six different 
symmetries analyzed (the corresponding values of the coupling constants are those reported above). In all the 
symmetries apart from the s-wave one, a conductance excess below the gap is obtained. However, in all cases 
the normalized conductance is zero at zero voltage. As we will discuss later, the zero-bias anomaly which is often 
experimentally observed could in fact be reproduced, in the framework of the Eliashbcrg theory, by taking into account 
the scattering from impurities, which has been disregarded here. 

As anticipated in the previous section, our choice of directly solving the real-axis EEs can be motivated by the 
fact that the analytical continuation of the solutions A(iw„) and Z(iuj n ) of the EEs in the imaginary-axis formulation 
is correct only at very low temperature. The meaning of this 'very low' can be clarified by comparing the results 
of the two approaches at various temperatures. The result of such a comparison is that, even at the liquid helium 



symmetry 


A is 
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s 


3.15 




95 




d 


2 


2.3 




95 


s+id 


2 


2.3 


66 


95 


s+d 


1.25 


1.74 


40 


95 


anis. s 


2 


2.29 


63 


95 


ext. s 


1 


1.54 


32 


95 



TABLE I: Values of the coupling constants used in the calculations of Section III and corresponding critical temper- 
atures for all the gap symmetries considered. 
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FIG. 6: Theoretical normalized con- 
ductances curves at T =4 K in various 
gap symmetries, for [i* = and in the 
presence of scattering from impurities 
in the non-unitary limit (r = 40 meV 
and C = 4). 



temperature (at which most of the experimental data are obtained), the an- 
alytical continuation gives rise to some sensible deviations from the real-axis 
solution. As shown in Fig. 2, at a lower temperature (T = 2 K), the agree- 
ment between the theoretical G n (V) curves obtained by the two procedures 
is good for the d-wave, s + id-wave and anisotropic s-wave symmetries. In 
the extended s case, the 'imaginary-axis' curve agrees very well with the 
'real-axis' one, except that in the low-voltage regime. In the (s + d)-wave 
symmetry, the peak of the imaginary-axis solution is broader than that of 
the real-axis one, and displaced of about +8 mV. Finally, in the s-wave case, 
the analytical continuation does not work well, since a shift of +3 mV of the 
conductance curve is produced. Further checks demonstrate that the agree- 
ment gets worse at the increase of A, while for sufficiently small values of the 
coupling constant (A < 2) the two curves coincide, as found in Ref. [13|. 

Let us now discuss what happens to the conductance curves when the 
temperature is increased. Figures 3, 4 and 5 show the normalized tunneling 
conductance curves calculated for all the symmetries at T = 2, 4, 40, 80 K. In 
general, as expected, increasing the temperature results in a smoothing and 
broadening of the peaks. The smoothing is particularly evident in the case 
of the anisotropic s symmetry, in which the low-voltage peak, which is well 
defined at T = 2 K, is much less evident at 4 K (its height is reduced from 
1.5 to 1.25). At the increase of the temperature the curves corresponding 
to mixed symmetries become more and more similar to one another and to 
the d-wave curve. This is due both to the thermal broadening and to the 
fact that, with our choices of Ai S and A an , T c ls is always less than T c an and 
therefore the isotropic gap component disappears before the anisotropic one. 
The s-wave curve remains instead clearly distinct. 
The theoretical curves presented and discussed so far were obtained without inserting in the EEs the additional 
term taking into account the scattering from non-magnetic impurities. When this effect is allowed for in solving the 
real-axis EEs, one finds that in the general non-unitary case (C 7^ 0) it mostly affects the d-wave component of the 
order parameter. By comparing the curves shown in Fig. 6 to those of Fig. 1, we can see that the peak of the d-wave 
tunneling conductance curve is lowered and broadened, while the s + id curve becomes practically indistinguishable 
from an s-wave one. In the s + d and extended s cases the low-energy peak disappears and a zero-bias is produced, 
which is greater in the latter case. Finally, in the anisotropic s case the two peaks are closer than in the absence of 
impurities and the gap is shifted to the left, but the general shape above and below the gap is conserved. 

In the unitary limit (C <C 1) the presence of impurities has very different effects in the various symmetries, as 
shown in Figs. 7, 8 and 9 that report the curves obtained with C = 0.1 and T = 1, 5, 10 meV. Also in this case the 
results indicate that, as expected, the main effect of non-magnetic impurities is to suppress the d-wave contribution. 
In the d-wave case, a small amount of impurities makes the peaks of the tunneling conductance curves to be lowered, 
broadened and shifted to the left, and for T > 10 meV the superconductivity disappears. In the same limit, the s+id 
and the anisotropic s curves become very similar to a pure s-wave, but the peak is shifted toward lower voltages. 
Finally, the peak of the extended s and s+d curves is lowered and broadened and the superconductivity is easily 
suppressed. 

An interesting feature of the theoretical G n (V) curves that occurs for suitable values of the parameters is a sup- 
pression of the tunneling conductance at about twice the gap value (usually referred to as the 'dip') and a consequent 
enhancement of conductance which occurs at about three times the gap value ('hump'). This feature, which has 
probably its origin in the non-linear character of the EEs, can be theoretically obtained at least in two ways. In 
a previous paper |l6| ] we showed that in s- and d-wave symmetries it can result from the finiteness of the energy 
band. In Ref. ]i~7[ , instead, 'dip' and 'hump' are obtained in the s-wave case by using a very large electron-boson 
coupling constant. Here we extend this result to the d-wave symmetry. Actually, in the /1* = case large values of 
the coupling constants would increase the critical temperature well beyond its physical value. To conciliate a very 
strong coupling with the correct T c , the maximum boson energy must be reduced, that is, the electron-boson spectral 
function must be confined in the low-energy region and, consequently, LOi og must be kept sufficiently small. Greater 
(and more plausible) values of o>i og can be obtained if one also takes fi* =/= 0. 

In Fig. 10 we report the theoretical d-wave conductance curves obtained with [i* = 0.2 and ui c — 500 meV for 
increasing values of the isotropic coupling constant: Aj S = 3,4,5. The corresponding anisotropic coupling constants 
(A an = 3.78,4.66,5.54) have been chosen so as to have T c calc = 95 K. The electron-boson spectral function has the 
form a 2 F(Q) is = bfl 2 exp(-VL / c) where c =8 meV and b =0.0234, 0.0313, 0.0391 meV~ 2 for A is = 3,4,5. In all cases 
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FIG. 7: Theoretical normalized con- 
ductance curves for fi* = in the d- 
wave symmetries at T= 4 K with im- 
purities in the unitary limit. 



FIG. 8: Theoretical normalized con- 
ductance curves for /i* = in the s+id 
and s + d symmetries at T= 4 K with 
impurities in the unitary limit. 



FIG. 9: Theoretical normalized con- 
ductance curves for fj* — in the 
anisotropic s and extended s symme- 
tries at T= 4 K with impurities in the 
unitary limit. 



c^iog = 12.2 meV, which is consistent with the energy of the first peak in the phonon spectral function of HTSC. It is 
clearly seen that the conductance curves reported in the figure show a 'dip' and a 'hump' at approximately 2A pea k/e 
and 3A P eak/e. 

Features very similar to those discussed so far have been indeed recently observed in the experimental conductance 
curves of some HTSC and reported in literature |Q. Notice that some possible explanations have already been 
proposed for these effects |2(| . 



IV. AN ILLUSTRATIVE APPLICATION 



Up to now we have dealt only with theoretical predictions, given by the Migdal-Eliashberg theory, concerning the 
conductance curves corresponding to different gap symmetries. We have discussed the effect of the temperature and 
of the scattering by impurities. 




25 50 75 100 125 150 175 200 

Voltage (mV) 



FIG. 10: Theoretical d-wave conductance curves calculated at T — 4 K for increasing values of the coupling constants. A 
suppression of the conductivity at V = 2A pca k/e ('dip') and an enhancement at l/ = 3A poa k/e ('hump') appear due to the strong 
electron-boson coupling. The inset shows the isotropic electron-boson spectral function a 2 F(Q)i B in the A s = 3 case. In the 
other cases the curve is simply scaled along the vertical axis. The anisotropic spectral function is a 2 F(Q) a , n = g ■ a 2 F(Q)i s . 
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Now we can give an example of how these calculations can be applied to a real physical problem. Actually, 
the Eliashberg theory has been successfully used to describe most of the low-T c superconductors by using a simple 
s-wave gap symmetry, thus it is not worthwhile to use these compounds to test our model. Moreover, we have 
performed our calculations in the 2-dimensional case, and therefore we should refer to layered materials. We can 
attempt a comparison of our theoretical predictions to tunneling conductance curves of SIN junctions involving high- 
T c superconductors (HTSC). Of course in applying this theory to HTSC, we must take into account that: i) a 2 n F(£l) 
doesn't necessarily have to be multiple of af s F(Q,); ii) fi* is probably different from zero; hi) as well known, the Fermi 
surface (FS) of these materials is not a circle; iv) the bandwidth is finite; v) the conductance of the junctions generally 
depends on geometrical factors which have been disregarded here. Actually, both a 2 F(fl, <f>, <j>') and ji* are unknown, 
and our assumptions are the simplest possible. The hypothesis that the FS is a circle is as well an oversimplification, 
but it allows obtaining a more handy model. As far as the bandwidth is concerned, we can presume on the basis of 
recent results that in optimally-doped HTSC e-p > 400 meV and therefore the effects of the finiteness of the bandwidth 
are not so important |16| ]. 

In general, the very complex phenomenology of HTSC yields to the conclusion that many features of these com- 
pounds are probably beyond the limits of the Eliashberg theory. Then, we don't expect this theory to provide an 
explanation for high-T c superconductivity. We will simply try to use our model to obtain, with suitable choices of the 
parameters, theoretical curves in agreement with experimental tunneling data in optimally-doped Bi 2 Sr 2 CaCu 2 0s+,5 



(BSCCO) g§, YBa 2 Cu 3 7 -5 (YBCO) gj, Tl 2 Ba2Cu0 6 + 4 (TBCO) |22| and HgBa 2 Cu0 4+ 5 (HBCO) |§ single 
crystals. Of course, the spectral function crF(Q, (j>, 4>') is unknown. As already pointed out, in the case of BSCCO 
we can take, according to eq.(ff^): 

a 2 F(n) is = ^a 2 F(n) Bi2212 



L 511 and for the other materials 



a 2 F(Q) is = ^G(fi) (16) 



where G(fi) is the phonon spectral density determined by neutron scattering p4[ , A is the corresponding coupling 
constant and Ai S is a free scaling parameter which must be adjusted to fit the experimental data. According to our 
previous assumptions, the anisotropic component of the spectral function is given by a 2 F(fl) an —g ■ a 2 F(fl) ls , where 
the constant g is another free parameter. Since in most of the experimental curves a large zero-bias conductance is 
present, we will solve the EEs (@)-(R) including the term which describes the scattering from impurities. 

Before going on with the discussion, let us focus for a while on the experimental data reported for BSCCO. As can 
be observed in Fig. 10, the break-junction tunneling measurements reported here give a conductance peak at an energy 
Apeak = eVp ea k = 25 meV, in good agreement with other break-junction data appeared in literature p5[ . However, 
this value is much smaller than that measured by other techniques, namely the scanning tunneling spectroscopy, which 
gives Ap ea k = 37 -j- 43 meV |27], |2^]. The reason for these discrepancies is not clear. On the other hand, we have 
already shown ^] that a large gap value such that measured by STM cannot be reconciled, in the framework of the 
Eliashberg theory, with the observed T c . On the contrary, the gap obtained from our data can be well reproduced 
by our model. For this reason, in the following we will choose to analyze our break-junction experimental data on 
BSCCO. 

We can now discuss the results of the comparison of experimental conductance curves to the theoretical ones 
calculated in the d and s+id symmetries. The choice of these symmetries is motivated by many experimental evidences 
for the existence of a d-wetve component of the order parameter in most of the superconducting copper-oxides. The 
agreement of theoretical and experimental curves is surprisingly good despite all the limitation cited above and the 
roughness of the model. In particular, the theoretical curves which give best results are those obtained in the s+id 
symmetry in the cases of TBCO and HBCO, and in the <i-wave symmetry in the other cases. Incidentally, it is 
worthwhile to stress that a good agreement can be obtained for different values of Ai s and g and, therefore, what we 
present here is only one of the possible choices. 

The results of the comparison are shown in Fig. 11, and the values of the parameters are reported in Table II. 

In the cases of BSCCO and YBCO the theoretical curves follow well the behaviour of experimental data, and also 
the calculated T c agrees with the measured one. In both these cases, due to the low zero-bias anomaly (~ 0.2), the best 
result is obtained with a very small impurity content (r = 0.5 meV and Y — 0.6 meV, respectively). The zero-bias 
has approximately the same value also in the cases of TBCO (0.19) and HBCO (0.23). However, in these cases the 
best agreement is obtained with Y = 0, because in the presence of impurities the peak of the theoretical curve results 
too low with respect to the experimental one. It is easily seen that the agreement is better for TBCO than for HBCO. 
In this latter case, to obtain an acceptable result we were forced to use parameters that give a calculated critical 
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FIG. 11: Comparison of the theoretical normalized tunneling conductance curves to the experimental data obtained on 
BSCCO (a) @, YBCO (b) TBCO (c) Q, HBCO (d) The symmetry of the order parameter used in the calculations 
is d-wave in cases (a) and (b), and s+id-wave in cases (c) and (d). For details on the values of the parameters see Table II. 



temperature T c calc = 143 K which is much greater than the measured one (T c exp = 97 K), while in all the other cases 

^calc ^cxp 

Let us conclude with a brief discussion of the assumption //* = we made at the beginning of this paper. This 
simplifying assumption could appear poorly adequate to describe the HTSC, where the Coulomb pseudopotential is, 
as well known, different from zero. Actually, the main effect of each component of (i* (isotropic and anisotropic) is 
to change the scale of the corresponding coupling strength. Then, almost the same tunneling curves can be obtained 
with fi* = and given values of the coupling constants, let's say and Ain > or with /i* s and /i* n different from zero 
and coupling constants A^ > A;*' and A^ > Ain ■ 



V. CONCLUSIONS 

We calculated the theoretical normalized tunneling conductance of SIN junctions for six different symmetries of the 
superconducting order parameter, by solving the real-axis EEs both in the presence and in the absence of scattering 
from non-magnetic impurities. We demonstrated that, in spite of being a much more complicated and time-consuming 
procedure, this approach is preferable to the analytical continuation of the solution of the imaginary-axis EEs, espe- 
cially if a comparison of the theoretical predictions with good experimental data is required. 





BSCCO 


YBCO 


TBCO 


HBCO 




2 


3.2 


2 


3 


A an 


2.32 


3.12 


2.2 


3.6 


r (meV) 


0.5 


0.6 








c 


0.1 


0.1 






T c cxp (K) 


93 


90 


91 


97 


r calc (jg 


93 


90 


91 


143 



TABLE II: Values of the parameters used in the calculations to obtain the theoretical normalized conductance curves 
that are compared to the experimental data for different high-T c superconductors. 
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For each of the symmetries considered, we also studied the temperature dependence of the tunneling conductance 
curves, showing that at the increase of temperature the information about the pairing symmetry is progressively lost. 

Finally, we discussed an example of how our calculations can be applied to a physical system, by comparing the 
theoretical curves to experimental tunneling data obtained on BSCCO, YBCO, TBCO and HBCO. Of course, it must 
be borne in mind that the applicability of the Eliashberg theory to these materials is controversial, and the reliability 
of the results is limited by the very complex phenomenology of high-T c superconductors. 

Unexpectedly, the results of such a comparison are fairly good. The theoretical curves reproduce some of the unusual 
properties of the tunneling curves of HTSC: the high T c , the broadening of the conductance peak, the conductance 
excess below the gap, the zero bias and also the presence of a 'dip' and a 'hump' at 2A/e and 3A/e, respectively. The 
results also seem to confirm, for these materials, the plausibility of a d-wave gap symmetry, or at least of a dominant 
ci-wave component of the order parameter. 

Even though alternative approaches to high-Tc superconductivity are nowadays preferred, maybe all these results 
indicate that, as far as the tunneling conductance is concerned, the Eliashberg theory mimics a (still unknown) theory 
of high-Tc superconductivity, whose mechanism could at least partially differ from the electron-boson coupling. 
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